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Abstract 

A  steady-state  mathematical  model  for  the  ordered  cathode  of  proton  exchange  membrane  fuel  cells  is  developed  to  investigate  the  dependence  of 
the  cathode  performance  on  the  structural  parameters  of  the  catalyst  layer.  The  model  is  based  on  the  governing  equations  for  oxygen  concentration 
and  potentials  of  the  membrane  and  the  solid  phase,  coupled  by  Tafel  relation  for  the  oxygen  reduction  reaction  kinetics.  The  cathode  current 
density  optimization  at  a  given  electrode  potential  is  presented  with  respect  to  nano-thread  radius,  porosity,  platinum  mass  percentage,  thickness, 
Nafion  volume  fraction  and  platinum  loading  of  the  catalyst  layer.  The  simulation  results  suggest  that  small  nano-thread  radius  is  preferred.  Except 
for  quite  low  values  as  well  as  thin  catalyst  layers,  porosity  and  platinum  mass  percentage  have  minor  effects  on  cathode  optimization.  The  cathode 
performance  depends  strongly  on  the  catalyst  layer  thickness  and  additional  attention  should  be  paid  to  a  thinner  catalyst  layer.  The  cathode  can  be 
efficiently  optimized  by  increasing  the  highly  sensitive  parameters,  Nafion  volume  fraction  and  platinum  loading,  to  a  suitable  value  which  must 
avoid  significant  loss  of  oxygen  transport. 
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1.  Introduction 

The  performance  of  proton  exchange  membrane  fuel  cells 
(PEMFCs)  is  known  to  depend  strongly  on  cathode  catalyst 
layers  due  to  the  usually  poor  kinetics  of  the  O2  reduction  reac¬ 
tion  (ORR)  and  transport  limitations  of  reactants.  Also,  the  high 
platinum  loading  needed  to  maintain  sufficient  reaction  rates  in 
cathodes  is  one  of  the  cost  factors  hampering  commercialization 
of  PEMFCs.  Thus,  in  order  to  improve  the  PEMFC  performance 
and  the  catalyst  utilization,  it  is  essential  to  analyze  and  optimize 
the  structure  and  performance  of  the  cathode  catalyst  layer. 

For  efficient  functioning,  at  least  four  ingredients  are  indis¬ 
pensable  in  a  catalyst  layer  (CF):  an  electron  conductor  (typi¬ 
cally  some  kind  of  carbon),  a  proton  conductor  (usually  Nafion 
polymer),  platinum  catalysts  and  gas  pores.  In  a  typical  catalyst 
layer  produced  by  a  number  of  conventional  techniques  [1-3], 
the  four  ingredients  are  randomly  distributed  which  results  in 
significant  amount  of  the  wasted  platinum  and  low  electrochem¬ 
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ical  performance  of  the  catalyst  layer  [4,5].  Theoretically,  an 
idealized  catalyst  layer  would  be  of  the  property  that  all  platinum 
particles  would  be  easily  accessible  for  hydrogen  or  oxygen, 
proton  and  electron  at  the  same  time  [6].  Recently,  Middel- 
man  [7]  suggested  a  novel  ordered  catalyst  layer  (OCF)  that  can 
achieve  partially  the  characteristic  of  an  idealized  catalyst  layer. 
The  OCF  features  a  highly  oriented  morphology  to  enhance 
catalyst  utilization  and  performance.  To  create  the  ordered  struc¬ 
ture,  the  mobility  of  the  catalyst  layer  was  increased  with  high 
temperatures  and  chemical  additives.  Then  an  electric  field  was 
employed  as  the  driving  force  to  orient  the  strands.  Middelman 
suggests  that  this  method  could  increase  the  Pt  utilization  to 
almost  100%  in  addition  to  the  20%  increase  in  voltage. 

Thorough  understanding  of  the  relation  between  structural 
characteristics  and  performance  is  important  to  optimize  the 
OCF  further.  Numerical  investigation  is  inexpensive,  less  time 
consuming  compared  to  experimental  methods,  in  particular, 
when  systematic  investigation  involving  more  than  one  parame¬ 
ter  is  conducted.  The  goal  of  this  study  is  to  analyze  and  optimize 
the  ordered  cathode  performance  with  respective  to  some  impor¬ 
tant  OCF  design  parameters:  the  nano-thread  radius,  porosity, 
the  mass  percentage  of  Pt  catalyst  supported  on  the  carbon  black, 
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thickness,  the  Nafion  volume  fraction,  and  the  platinum  loading 
by  a  numerical  method.  For  this  purpose,  a  numerical  model 
of  the  ordered  cathode  in  PEMFCs  is  presented  by  including 
the  processes  of  the  electrochemical  kinetics,  the  oxygen  trans¬ 
port,  and  proton  and  electron  transfers.  The  results  highlight 
the  interdependence  of  the  cathode  performance  and  the  OCL 
design  parameters.  It  is  hoped  that  these  results  will  be  helpful 
to  the  optimization  and  design  of  the  OCL. 


the  proton  transfer  in  the  CL  and  the  PEM;  the  electrochemical 
reaction  of  oxygen  reduction  on  the  catalyst  surface. 

Oxygen  transport  in  the  GDL  and  the  CL  is  mainly  due  to  the 
oxygen  concentration  gradient  and  approximately  described  by 
Fick’s  law  with  the  following  equation  forms  [8]: 

For  GE>L,  A(Dg-i^)=0  0) 


2.  Model  description 


For  CL, 
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dx 


(  T\ CL  dc02,g^\ 
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The  ordered  cathode  model  presented  here  is  a  one¬ 
dimensional,  isothermal  and  steady-state  model.  The  catalyst 
layer  is  assumed  to  be  a  homogeneous  medium  consisting  of 
oriented  carbon  nano- threads  which  are  covered  by  Nafion  elec¬ 
trolyte  films  and  then  separated  by  parallel  gas  micropores.  The 
catalyst  is  assumed  to  be  uniformly  distributed  on  the  carbon 
nano-thread  surface.  The  Nafion  electrolyte  is  assumed  to  be 
fully  humidified  and  have  constant  proton  conductivity.  The 
computation  domain  of  the  ordered  cathode,  as  shown  in  Fig.  1, 
is  composed  of  the  proton  exchange  membrane  (PEM),  the  OCL 
and  the  gas  diffusion  layer  (GDL).  The  spatial  coordinate  v  is 
defined  so  that  the  positive  direction  points  from  the  membrane 
to  the  electrode  with  its  origin  located  at  the  membrane/anode 
interface. 

2.7.  Governing  equations 

The  interconnected  physical  and  chemical  processes  mod¬ 
eled  to  describe  the  behavior  of  the  cathode  are  as  follows:  the 
oxygen  diffusion  in  the  GDL  and  the  CL;  the  oxygen  dissolution 
at  the  gas-electrolyte  film  interface;  the  diffusion  of  dissolved 
oxygen  in  the  Nafion  electrolyte  film  covered  on  the  catalyst;  the 
electron  conduction  in  the  cathode,  both  in  the  GDL  and  the  CL; 


where  co2,g  is  the  oxygen  concentration  in  the  gas  phase, 

and  D%\  eff  the  effective  diffusion  coefficients  of  oxygen  in  the 
respective  GDL  and  CL,  i  the  volumetric  current  density,  n  the 
number  of  electrons  transferred  in  the  ORR,  and  F  the  Faraday’s 
constant. 

To  account  for  gas  diffusion  in  the  porous  GDL,  the  oxygen 
diffusion  coefficient  has  been  modified  to  the  effective  diffusiv- 
ity  through  the  commonly  employed  Bruggemann  relation  [9]: 

^O^eff  =  ^02,g£GDL  (3) 

where  £gdl  is  the  porosity  of  the  GDL,  Z)o2,g  the  bulk  oxygen 
diffusion  coefficient  in  the  gas  phase.  In  view  of  the  parallel  gas 
pores,  the  oxygen  diffusion  coefficient  in  the  CL  need  not  be 
corrected  by  tortuosity  factor  and  is  calculated  as  follows: 

^O^eff  =  A}2,g£CL  (4) 

where  £cl  is  the  porosity  of  the  CL  and  can  be  determined  by: 

£CL  =  1  -  £n  ~  £s  (5) 

where  £n  is  the  Nafion  volume  fraction,  and  £s  is  the  volume 
fraction  of  the  solid  phase  (catalyst  plus  carbon)  in  the  CL  and 


Fig.  1.  Schematic  diagram  of  the  PEMFC  cathode  model  and  a  typical  cross-section  of  the  OCL. 
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is  given  by  [10]: 


£s 


1  1  -  mpPt  \  mlpt 

PPt  mpptPc  )  LCl 


where  ppt  and  pc  are  the  mass  densities  of  the  Pt  and  carbon 
black,  respectively,  mppt  the  mass  percentage  of  Pt  catalyst  sup¬ 
ported  on  the  carbon  black,  and  mlpt  is  the  mass  loading  of  Pt 
per  unit  area  of  the  CL. 

The  dissolved  oxygen  concentration  at  the  gas-electrolyte 
interface  is  assumed  to  be  in  equilibrium  with  its  concentration 
in  the  gas  phase  and  is  described  by  Henry’s  Law  [11]: 


c02,m 


co2,g 

Ho2 


where  co2,m  is  the  dissolved  oxygen  concentration  in  the  Nation 
electrolyte  membrane  and  Ho2  the  Henry’s  constant  calculated 
from  the  empirical  correlations  [12]: 

1  /  666  \ 

Hq2  =  —  exp  ( — —  +  14.1  )  (8) 


where  T  is  the  cathode  temperature,  R  the  gas  constant. 

The  dissolved  oxygen  must  diffuse  cylindrically  from  the 
gas-electrolyte  interface  to  the  catalyst  sites  for  ORR.  This 
process  is  governed  by  concentration  gradient  and  yields  the 
following  equation  [13]: 
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where  Ra  and  R\>  are  the  inner  and  outer  radii  of  the  cylindri¬ 
cal  electrolyte,  A  the  catalyst  specific  area,  and  k  is  defined  as 
follows: 


(13) 


Since  the  nano-thread  is  cylindrical,  the  outer  radius  of  the  cylin¬ 
drical  electrolyte  Rb  is  related  to  the  inner  radius  Ra  by 


£N  +  £S 
£S 


The  specific  active  surface  area  A  is  calculated  by 


(14) 


(15) 


where  A0  is  the  catalyst  surface  area  per  unit  mass  of  the  catalyst. 

The  governing  equations  for  the  proton  migration  in  the  mem¬ 
brane  and  the  CL  can  be  derived  from  the  continuity  equation 
and  Ohm’s  law  as  follows  [17]: 


in  the  membrane, 


(16) 


where  Do2,m  is  the  oxygen  diffusion  coefficient  in  the  electrolyte 
membrane,  r  the  radial  coordinate  in  the  cylindrical  electrolyte. 

The  rate  of  the  ORR,  which  is  a  function  of  the  oxygen 
concentration  and  the  local  overpotential,  is  described  by  the 
phenomenological  Butler- Volmer  theory.  In  view  of  the  typical 
current  density  in  operating  PEMFCs,  only  the  forward  current 
term  in  the  Butler- Volmer  expression  is  adopted  [14]: 

*s  =  i°  (~^r)  exP  yoc))  (io) 

where  is  is  the  rate  of  the  ORR,  fi  the  exchange  current  den¬ 
sity,  cref  the  reference  oxygen  concentration,  aY  the  cathodic 
transfer  coefficient,  Voc  the  thermodynamic  equilibrium  poten¬ 
tial,  and  cps  and  cpm  are  the  potentials  of  the  carbon  phase  and 
the  membrane,  respectively.  The  exchange  current  density  can 
be  correlated  to  the  experimental  data  with  temperature  in  the 
following  expression  [15]: 

n  4001 

logjoO'0)  =  3.507  -  —  (11) 

The  equations  of  the  oxygen  dissolution  at  the  gas-electrolyte 
interface,  the  cylindrical  diffusion  of  oxygen  in  the  Nation  film 
and  the  ORR  rate  (Eqs.  (7),  (9)  and  (10))  can  give  the  overall 
volumetric  current  density  /,  which  is  expressed  in  terms  of  the 
oxygen  concentration  at  the  outer  surface  of  the  cylindrical  nano¬ 
thread,  the  membrane  potential  and  the  solid  phase  potential 


intheCL,  ±  ^ )  =  -i  (17) 

where  crm  is  the  bulk  proton  conductivity  of  the  Nation  mem¬ 
brane,  and  a^Leff  is  the  effective  proton  conductivity  in  the  CL. 

The  electron  transfer  in  the  GDL  and  the  CL  can  also  be 
derived  from  the  mass  balance  and  Ohm’s  law  and  described  by 
the  following  governing  equations  [17]: 

in  the  GDL,  ±  (<®  £)  =  0  (18) 

'"lheCL-  <19) 

where  erf  and  are  the  effective  electronic  conductivities 
in  the  CL  and  the  GDL. 

Similar  to  the  case  of  the  oxygen  diffusion  coefficient,  the 
effective  proton  and  electron  conductivities  in  the  CL  and  the 
GDL  can  be  determined  in  the  following  forms: 

in  the  GDL,  (20) 

intheCL,  =  orseCL,  ^ni'eir  =  °m£CL  (21) 

2.2.  Boundary  conditions 

The  model  Eqs.  (1)  (2)  and  (14)— (17)  are  solved  subject  to 
the  following  boundary  conditions: 
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(a)  Boundary  conditions  for  the  oxygen  transport: 


atx  -  Lm,  2,g  -  0; 

ax 

(22) 

atx  =  Lm  +  Lql  +  ^gdl,  <A>2,g  =  co2,g; 

(23) 

(b)  Boundary  conditions  for  the  proton  transfer: 

atx  =  0,  (Pm  =  0; 

(24) 

dcpm 

atx  =  Lm  +  LCl,  ,  =0; 

ax 

(25) 

(c)  Boundary  conditions  for  the  electron  transfer: 

atx  =  Lm  +  Lql  +  Tgdl,  <Ps  = 

(26) 

dcp, 

atx  =  Lm.  .  =0; 

ax 

(27) 

3.  Results  and  discussion 

The  model  equations  together  with  the  boundary  conditions 
listed  in  Section  2  are  solved  with  the  finite  element  method 
and  the  program  is  written  with  the  C  language  by  modifying 
the  AFEPack  code  developed  by  Peking  University.  Each  gov¬ 
erning  equation  is  solved  throughout  the  entire  domain,  which 
is  accomplished  by  setting  zero  equation  coefficients  in  phys¬ 
ically  invalid  subdomains.  Sensitivity  analysis  shows  that  316 
elements  in  the  v-direction  provide  adequate  resolution  and  the 
same  element  number  is  used  for  all  simulations. 

The  operating  parameters  and  physical  properties  for  the 
base-case  simulation  are  listed  in  Table  1.  The  void  fraction 
scl ,  the  Nation  volume  fraction  £n,  the  platinum  loading  mlpt, 
the  mass  percentage  of  Pt  catalyst  mppt,  the  catalyst  layer  thick¬ 
ness  Lcl,  and  the  nano-thread  radii  La  and  R b  are  related  by  Eqs. 
(5),  (6)  and  (14).  Consequently,  for  a  given  type  of  catalyst  layer, 
only  five  (must  include  La  and/or  Lb)  out  of  the  seven  parameters 
(£cl>  £n,  mlpt,  mppt,  La,  R b  and  Lql)  need  to  be  specified,  and 
the  sixth  and  the  seventh  as  well  as  £s  can  be  determined  from 
Eqs.  (5),  (6)  and  (14).  In  the  following  parametric  studies,  £ n, 
mlpt,  mppt,  Lql  and  therefore  £cl  are  fixed  for  the  case  of  the 
nano-thread  radius;  except  that,  four  out  of  the  seven  parameters 
are  fixed  and  the  other  two  vary  with  the  investigated  parameter, 
satisfying  the  constraint  conditions  Eqs.  (5),  (6)  and  (14). 

3.1.  Effect  of  the  nano-thread  radius 

Fig.  2  shows  changes  of  the  current  density  with  the  nano¬ 
thread  radius  at  a  typical  electrode  potential.  It  is  observed 
that  the  current  density  decreases  by  increasing  the  nano-thread 
radius,  with  a  little  more  pronounced  effect  for  higher  nano¬ 
thread  radii.  Moreover,  variations  of  the  current  density  are 
almost  identical  for  different  CL  thicknesses.  As  the  nano-thread 
radius  increases  from  10  to  190  nm,  the  total  decrease  of  the  cur¬ 
rent  density  is  about  100mA cm-2.  This  suggests  that  thinner 
nano-threads  are  preferred  for  the  efficient  operation  of  the  OCL. 
However,  in  the  design  of  the  OCL,  the  nano-thread  radius  has 


Table  1 


Operating  and  structural  parameters  used  in  the  base-case  simulation  [18,19] 


Parameter 

Symbol 

Value 

Surface  area  per  unit  mass  of  the  catalyst 

A° 

112 

(m2g  l) 

Oxygen  concentration  at  the  GDL  surface 

c02,g 

7.46 

(molm  3) 

Reference  oxygen  concentration  (mol  m  3) 

1.2 

02  diffusion  coefficient  in  gas  phase 

£>o2,g 

2.837  x  10“5 

(m2  s  1 ) 

02  diffusion  coefficient  in  the  Nafion 

Do2,m 

6.0  x  10“10 

(m2  s  1 ) 

Faraday’s  constant  (C  mol-1) 

F 

96485 

The  CL  thickness  (m) 

CcL 

1.5  x  10"5 

The  GDL  thickness  (m) 

Tgdl 

2.5  x  10“4 

The  PEM  thickness  (m) 

Lm 

1.0  x  10“4 

Mass  percentage  of  platinum  supported  on 

mppt 

0.2 

carbon 

Catalyst  mass  loading  per  unit  area  of  the 

mlpt 

0.3 

cathode  (gm  2) 

Number  of  electrons  transferred  in  the 

n 

4 

reaction 

Gas  constant  (J  K-1  mol-1) 

R 

8.314 

The  carbon  nano-thread  radius  (m) 

Ra 

3.9  x  10"8 

Temperature  (K) 

T 

343 

Thermodynamic  equilibrium  potential  (V) 

Foe 

1.19 

Transfer  coefficient  of  the  ORR 

ar 

0.58 

The  GDL  porosity 

£gdl 

0.4 

The  Nafion  volume  fraction 

£n 

0.25 

Mass  densities  of  the  platinum  (gm-3) 

pPt 

21.5  x  106 

Mass  densities  of  the  carbon  black  (gm-3) 

PC 

2  x  106 

Bulk  electron  conductivity  in  the  CL  (S  m-1) 

72700 

Bulk  electron  conductivity  in  the  CL  (S  m-1) 

7 

to  be  carefully  selected  because  too  thin  nano- threads,  in  prac¬ 
tice,  may  result  in  the  reduction  of  the  platinum  loading  and  thus 
lower  the  CL  performance.  Effects  of  the  nano-thread  radius  on 
the  current  density  can  be  ascribed  to  the  thickness  variation  of 
the  Nation  film  covered  on  the  catalyst  surface  (also  shown  in 
Fig.  2).  From  Fig.  2  and  also  Eqs.  (6)  and  (14),  the  Nation  film 
thickness  covered  on  the  catalyst  ( R b  —  La),  increases  linearly 
with  increasing  the  nano-thread  radius  when  £n,  mlpt,  mppt  and 
Lql  are  kept  constant,  which  causes  increases  in  oxygen  diffu- 


Fig.  2.  Effects  of  the  nano-thread  radius  on  the  current  density  at  0.6  V  and  the 
Nation  film  thickness  on  the  catalyst  surface. 
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Fig.  3.  Effects  of  the  CL  porosity  on  the  current  density  at  0.6  V  and  the  Nation 
film  thickness  covered  on  the  catalyst  surface. 


sion  resistance  through  the  electrolyte  film  and  hence  lowers  the 
current  density  at  a  given  electrode  potential. 


3.2.  Effect  of  the  porosity 

Fig.  3  shows  changes  of  the  current  density  with  the  porosity 
of  the  OCL  for  fixed  £n,  mlpt  and  Lcl-  The  current  density 
remains  almost  constant  except  for  the  case  of  very  low  porosi¬ 
ties  when  the  effective  oxygen  diffusion  coefficient  in  the  OCL 
is  lowest  according  to  Eq.  (4).  The  low  sensitive  performance 
dependence  on  the  porosity  can  be  attributed  to  the  improvement 
of  oxygen  transport  by  orientation  of  the  OCL.  Note  that  the  cur¬ 
rent  density  decreases  also  a  little  for  very  high  porosities.  This 
may  be  due  to  the  quick  increase  of  the  Nation  film  thickness 
covered  on  the  catalyst  (also  shown  in  Fig.  3)  for  a  fixed  Nation 
volume  fraction  in  the  OCL,  which  is  derived  from  Eq.  (5),  (6) 
and  (14): 


Rb  -  Ra  =  ( J  1  £CL  -  l)  R* 

\  v  i  —  £cl  ~  £n  J 

3.3.  Effect  of  the  platinum  mass  percentage 


(28) 


Fig.  4  shows  the  relationship  between  the  current  density  and 
the  mass  percentage  of  the  platinum  supported  on  the  carbon 


mpP t.  % 

Fig.  4.  The  current  density  with  the  mass  percentage  of  platinum  for  different 
CL  thicknesses  at  a  given  electrode  potential  (0.6  V). 


black  for  fixed  7?a,  £n,  mlpt  and  Lcl  -  In  the  simulated  range  of  the 
platinum  mass  percentage,  the  current  density  exhibits  distinct 
behaviors  for  different  CL  thicknesses.  For  a  thicker  CL,  the 
current  density  remains  almost  constant  as  the  platinum  mass 
percentage  increases;  while  for  a  thinner  CL,  the  current  density 
increases  sharply  and  then  decreases  slightly  with  increasing  the 
platinum  mass  percentage.  These  observations  can  be  explained 
by  changes  of  the  CL  porosity  and  the  Nation  film  thickness  as 
a  result  of  varying  the  platinum  mass  percentage. 

Fig.  5  gives  the  CL  porosity  as  a  function  of  the  platinum 
mass  percentage  and  the  platinum  loading  for  fixed  Ra,  £n  and 
Lcl-  Obviously,  as  the  platinum  mass  percentage  increases,  the 
CL  porosity  increases  slightly  except  near  the  lowest  platinum 
mass  percentage.  It  is  also  shown  in  Fig.  5  that  the  CL  porosity  is 
more  sensitive  to  the  platinum  mass  percentage  for  a  thinner  CL 
than  a  thicker  one.  These  results  agree  well  with  those  in  Fig.  4 
except  the  slight  decrease  of  the  current  density  at  high  platinum 
mass  percentages.  This  deviation  is  attributed  to  the  thickness 
increase  of  the  Nation  film  covered  on  the  catalyst  in  the  case  of 
fixed  Ra ,  £n  and  Lcl,  which  is  derived  from  Eqs.  (6)  and  (14)  and 
shown  in  Fig.  6.  In  this  figure,  the  thickness  of  the  Nation  film 
increases  with  increasing  the  platinum  mass  percentage,  which, 
in  turn,  increases  the  oxygen  transport  resistance. 
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Fig.  5.  Dependence  of  the  CL  porosity  on  the  platinum  mass  percentage  and  the  platinum  loading  for  two  CL  thicknesses. 
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Fig.  6.  Dependence  of  the  Nation  electrolyte  film  thickness  on  the  platinum  mass  percentage  and  the  platinum  loading  for  two  CL  thicknesses. 


3.4.  Effects  of  the  catalyst  layer  thickness 

Effects  of  the  CL  thickness  on  the  cathode  performance  for 
fixed  Ra,  £n,  mlpt  and  mppt  are  presented  in  Fig.  7.  It  is  seen 
that  there  is  an  optimal  CL  thickness  where  the  current  density 
can  reach  2.05  A  cm-2  for  the  0.2  platinum  mass  percentage. 
Furthermore,  the  current  density  is  much  more  sensitive  to  thin¬ 
ner  CLs  than  to  thicker  ones,  in  particular,  below  the  optimal 
thickness.  This  means  that,  in  the  CL  design,  thickness  needs 
careful  consideration  especially  for  thinner  CLs.  Effects  of  the 
CL  thickness  can  be  explained  by  variations  of  the  proton  and 
oxygen  transport  resistances.  As  the  CL  thickness  decreases,  the 
path  lengths  of  the  proton  and  oxygen  transport  are  shortened 
and  the  Nation  film  thickness  on  the  catalyst  shown  in  Fig.  8  is 
reduced,  which  improve  the  CL  performance.  However,  exces¬ 
sively  low  CL  thickness  results  in  extremely  low  CL  porosity 
(also  in  Fig.  8)  and  reduces  CL  performance. 

3.5.  Effects  of  the  Nafion  volume  fraction 

Effects  of  the  Nafion  volume  fraction  on  the  current  density 
are  presented  in  Fig.  9  for  different  CL  thicknesses  in  the  case 
of  fixed  Ra ,  mlpt  and  mppt,  and  in  Fig.  10  for  different  platinum 


Fig.  7.  The  current  density  at  0.6  V  electrode  potential  as  a  function  of  the  CL 
thickness  for  different  platinum  mass  percentages. 


Fig.  8.  The  CL  porosity  and  the  Nafion  film  thickness  on  the  catalyst  as  a  function 
of  the  CL  thickness. 

loadings  in  the  case  of  fixed  Ra,  mppt  and  Lcl-  From  the  two 
figures,  it  is  shown  that  as  the  Nafion  volume  fraction  increases 
the  current  density  increases  and  gradually  levels  off  to  a  high 
value.  Further  increases  of  the  Nafion  volume  fraction  result  in 
even  the  current  density  drop,  especially  for  the  cases  of  thinner 
CLs  and  high  catalyst  loadings.  These  results  can  be  explained 


Fig.  9.  Effects  of  the  Nafion  volume  fraction  on  the  current  density  at  0.6  V 
electrode  potential  for  different  CL  thicknesses. 
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Fig.  10.  Effects  of  the  Nation  volume  fraction  on  the  current  density  at  0.6  V 
electrode  potential  for  different  platinum  loadings. 

as  follows.  The  Nafion  volume  fraction  is  closely  related  to  the 
proton  transport  capability  of  the  CL.  When  the  Nafion  volume 
fraction  is  low,  the  poor  proton  transport  is  a  critical  factor  to  the 
CL  performance  and,  consequently,  increasing  the  Nafion  vol¬ 
ume  fraction  can  improve  the  current  density.  Simultaneously,  as 
the  Nafion  volume  fraction  increases,  the  Nafion  film  thickness 
on  the  catalyst  increases  and  the  CL  porosity  decreases,  which 
hamper  oxygen  transport.  These  two  factors,  especially  the  lat¬ 
ter,  will  result  in  a  fast  decline  of  the  current  density  for  very 
high  Nafion  volume  fractions.  Note  that  effects  of  the  Nafion 
volume  fraction  become  more  severe  for  thinner  CLs  and  higher 
catalyst  loadings,  which  is  obviously  attributed  to  the  higher  sen¬ 
sitivity  of  the  porosity  and  the  Nafion  film  thickness  on  thinner 
CLs  and  high  catalyst  loadings  as  shown  in  Fig.  5  and  Fig.  6, 
respectively.  It  should  also  be  noted  that  Nafion  is  assumed 
to  be  fully  humidified  in  the  present  simulation.  If  Nafion  is 
drier  and  hence  the  proton  conductivity  is  lower,  the  Nafion  vol¬ 
ume  fraction  will  have  more  pronounced  effects  on  the  OCL 
performance. 

3.6.  Effects  of  the  platinum  loading 

Fig.  1 1  shows  the  current  density  as  a  function  of  the  platinum 
loading  for  fixed  7?a,  £n,  mppt  and  Lcl-  Similar  to  the  case  of  the 
Nafion  volume  fraction,  the  current  density  increases  quickly  at 
the  beginning  and  then  drops  dramatically  as  the  catalyst  loading 
increases.  These  results  can  also  be  explained  by  the  compet¬ 
ing  effects  of  increasing  the  platinum  loading.  Higher  platinum 
loadings  give  larger  catalyst  active  surfaces  and  smaller  Nafion 
film  thicknesses,  which  increase  the  current  density  for  a  given 
electrode  potential.  For  too  high  platinum  loadings,  however, 
the  benefit  of  higher  catalyst  active  surfaces  does  not  make  up 
for  the  increased  oxygen  diffusion  resistance  as  a  result  of  the 
porosity  decrease,  which  results  in  the  performance  degradation 
of  the  CL.  Note  that  the  current  density  is  more  subtle  to  the 
platinum  loading  for  thinner  CLs.  This  is  also  attributed  to  the 
high  sensitivity  of  the  porosity  and  the  Nafion  film  thickness 
on  thinner  CLs.  From  these  results,  it  is  clear  that  not  too  high 


Fig.  1 1 .  Dependence  of  the  current  density  on  platinum  loading  for  different  CL 
thicknesses  at  0.6  V  electrode  potential. 


m/Pt/mg  cm'2 

Fig.  12.  Dependence  of  the  current  density  on  the  platinum  loading  for  different 
CL  porosities  at  0.6  V  electrode  potential. 

catalyst  loadings  are  preferred  especially  for  thinner  CLs  even 
regardless  of  the  high  catalyst  costs. 

Fig.  12  gives  effects  of  the  platinum  loading  on  the  current 
density  for  some  given  porosity  values  in  the  case  of  fixed  7?a, 
£n  and  Lcl-  With  increases  in  platinum  loading,  the  current 
density  is  quickly  enhanced  to  a  steady  value.  Furthermore,  the 
current  density  maintains  steady  even  at  1 .5  mg  cm-2  platinum 
loading  for  the  case  of  constant  porosity  whereas  the  current 
density  drops  quickly  at  0.5  mg  cm-2  platinum  loading  for  the 
case  of  variable  porosity  shown  in  Fig.  11.  This  demonstrates 
the  significance  of  maintaining  a  suitable  OCL  porosity. 

4.  Conclusions 

A  steady- state  numerical  model  is  developed  to  account 
for  the  ordered  cathode  in  PEMFCs  by  including  the  electro¬ 
chemical  reaction,  the  oxygen  transport  and  transfers  of  the 
two  charged  species — electron  and  proton.  Effects  of  the  cata¬ 
lyst  layer  structural  and  composition  parameters  on  the  ordered 
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cathode  performance  are  systematically  studied.  Based  on  the 

modeling  results,  the  following  can  be  concluded: 

(1)  The  nano-thread  radius  of  the  catalyst  layer  should  be  kept 
as  small  as  possible  to  minimize  the  oxygen  transport  resis¬ 
tance  through  the  Nation  film  covered  on  the  catalyst. 

(2)  Except  for  the  case  of  quite  low  values,  porosity  of  the  cat¬ 
alyst  layer  has  minor  effects  on  the  cathode  performance, 
which  results  from  the  oriented  gas  pores  enabling  fast  oxy¬ 
gen  diffusion  into  the  entire  depth  of  the  catalyst  layer. 

(3)  The  platinum  mass  percentage  of  the  catalyst  layer  has  visi¬ 
ble  effects  on  the  cathode  performance  just  when  the  catalyst 
layer  is  thin  and  the  platinum  mass  percentage  is  low.  Only 
in  this  case,  should  the  platinum  mass  percentage  be  care¬ 
fully  chosen  for  the  cathode  optimisation. 

(4)  The  cathode  performance  depends  largely  on  thicknesses  of 
the  catalyst  layer  and  additional  attention  should  be  paid  to 
a  thinner  catalyst  layer  during  optimization  in  order  to  avoid 
a  large  drop  in  current  density. 

(5)  The  cathode  performance  depends  highly  on  the  Nation  vol¬ 
ume  fraction.  The  optimal  performance  can  be  achieved  by 
adjusting  the  Nation  volume  fraction  to  a  suitable  high  value 
which  can  improve  proton  transfer  and,  at  the  same  time, 
maintain  a  large  enough  OCL  porosity  to  facilitate  oxygen 
transport  in  the  OCL. 

(6)  The  cathode  performance  is  much  more  sensitive  to  plat¬ 
inum  loading  than  to  the  other  parameters.  Increasing  the 
platinum  loading  can  increase  the  cathode  performance  effi¬ 
ciently.  However,  not  too  high  platinum  loading  is  preferred 
due  to  the  limitation  of  slow  gas  diffusion  in  the  case  of 
very  low  catalyst  layer  porosity,  especially  for  thinner  cata¬ 
lyst  layers. 
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